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We investigate the electronic Raman scattering in pure, quasi-one dimensional conductors with 
' density wave ground state. In particular, we develop the theory of light-scattering on spin and 

charge density waves, both conventional and unconventional. We calculate the electronic Raman 
response of the interacting electron system with a single, highly anisotropic conduction band. The 
OA ' calculation is carried out in the mean field approximation. Beside the quasiparticle contribution, the 

electron-electron interaction is also included on RPA level. The contribution of collective modes and 
q . the effect of Coulomb screening are investigated. In analogy with unconventional superconductivity, 

y-^ • the obtained Raman spectra - which are finite in the low temperature phase possessing a gap, and 

vanish identically in the normal state - show unique and strong dependence on the polarization 
, of the incoming and scattered light. We have found distinct, characteristic lineshapes, especially 

£N| ■ in the unconventional situation, depending on the various scattering geometries and the particular 

momentum dependence of the density wave order parameter. 

. PACS numbers: 75.30.Fv, 71.45.Lr, 78.30.-j, 72.15.Nj 
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^ . I. INTRODUCTION 

c3 . 

Electronic Raman scattering has been proven to be a valuable spectroscopic tool in identifying various low tem- 
perature phases of interacting electron systems. In early measurements on layered transition metal dichalcogenides 
like 2H-NbSe 2 , a low energy peak was attributed to the charge density wave amplitude mode pj , while at even lower 
temperatures the superconducting gap also showed up in the Raman spectra 0. A series of papers i nves tigated the 
Raman response of a system with competing spin density wave and superconducting instabilities [il l lol lop. Raman 
i_T experiments contributed significantly to the establishment of the d-wave nature of the order parameter in high tem- 
perature superconductors (HTSC) Q- It has also been applied recently in order to investigate the temperature and 
pressure dependence of the charge density wave amplitude mode in lT-TiSe2 . Superconduc ting a nd density wave 
condensates, both unconventional, are believed to be present in the underdoped cuprates lil lKIl ll . A theoretical 
analysis of this complex situation with respect to Raman scattering has also been attem pted [Hi . Recent work on 
Raman spectra in HTSC calls attention to the importance of density fluctuations as well |!3ll4ll5| . 
£NJ . The recent surge of in terest in u nconventional density waves (UDW) is mostly due to their potential applicability 
t— I ' in the pseudogap phase |!6ll7ll8| of HTSC materials. However pseudogap phases, and in general various kinds of 
hidden order, are detected in other substances as well, like in chalcogenides [l|| , in heavy fermion materials j2jj, and 
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in Bechgaard salts [21|. Since UDW's are natural candidates for explaining hidden order due to their momentum 
dependent gap structure they have been proposed to exist in URu9S i^> p 3l and in a-(ET)2 salts [24|. Recent 
calculations of magnetoresistance, thermoelectric power and Nernst effect |25| point out the possibility of UDW in 
(TMTSF) 2 PF 6 Hl| and CeCoIn 5 (27|. NMR results on Na . 7 CoC>2 are also consistent with the UDW scenario pjj . 

The aim of the present paper is to develop a theory of Raman scattering in pure quasi-one dimensional conductors 
with conventional, or unconventional density wave ground state. Basics of electronic Raman scattering and mean 
field treatment of density waves are given in Section II. The quasiparticle contribution to the light scattering intensity 
in various polarizations and gap structures are calculated in Section III. In Section IV we incorporate the effect of 
, electron-electron interaction. Namely we consider the collective contribution caused by the fluctuation of the order 
parameter in the usual RPA approximation, and also investigate the Coulomb screening. Finally Section V is devoted 
to our conclusions. 



II. ELECTRONIC RAMAN SCATTERING 

Light coupling to electrons via the vector potential A can be treated in second-order perturbation theory. The 
intensity of scattered light in a Raman experiment can be expressed j2t| as 

da n lu s 

___ = ro _S 77 (q, w ), (1) 

where r 2 , = e 2 /mc 2 is the Thomson radius, 0Ji,<\i and ui s ,q s are the energies and momenta of the incoming and 
scattered photon, respectively. Furthermore the energy and momentum transfer to the material are to — LUi — lu s and 
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q = q; — . The generalized structure factor iS 77 is related to the Raman response through the fluctuation-dissipation 
theorem 

Sry 7 (q,uj) = -[l + n(cj)]Imx T7 (q,w), (2) 
where n(oj) is the Bose function. The Raman response of the electron system measures "effective density" fluctuations 

Xrr(q»w)=i([p(q3,p(-q)])H/V; (3) 

where 

;5(q) = I]^ kC k+q, CT c k, CTI (4) 

k,(T 

V is the volume of the system, and since we are interested in the q — ► behavior of the x 77 susceptibility, we 
neglected the q dependence of the vertex 7k- Here c k(T , (ck, CT ) is the creation (annihilation) operator of an electron 
with momentum k and spin a in the single conduction band = — 2t a cos ak x — 2tb cos bk y — 2t c cos ck z with t a 3> tb,t c . 
Our system is based on an orthorombic lattice with lattice constants a, b,c towards the x,y and z directions. The 
strength of the scattering is determined by the momentum dependent function 7k called the Raman vertex, which 
has the form 



7k 



(c t c s ) I 1 y- / (k|pejfek)(fck|pejk) | (k|pe t |bk)(bk|pe s |k) \ 



where 6 stands for the band index of the electron excited out of the conduction band, and the corresponding states 
are |k&) and |k), respectively. In addition the polarization vectors of the incoming and scattered light are denoted by 
ei,e s . If the incoming and scattered light frequencies can be neglected in comparison to the optical band gap |30| . 
the Raman vertex is related to the inverse mass tensor 7 Q ^(k) = md 2 ey i /dk a dkp through the relation 

7 k = $>« 7ct/3 (k)ef , (6) 

a,/3 

which is widely known as the effective mass approximation. 

The retarded susceptibility of the effective density can be evaluated with analytical continuation from the Fourier 
transform of the corresponding r (imaginary time) ordered response Xj-y( c i: T ) — ~ (^r[p(q, r)p(— q)]) in the usual 
way. The one particle Green's function of the DW using Nambu's notation reads 



G(k,iu> n ) = - rdT(r r [*(k,r)vI/t(k)])e^ r , 
Jo 



(7) 



where the four component spinor field 

*(k,r) 



/ c k ,T(r) 

Ck-Q,tW 

Ck,i(r) 

Vck-Q,l( r )/ 



(8) 



is introduced to simultaneously cover the spin space and to treat the left- and right-moving electrons in momentum 
space in a convenient way. Q = (2kp, n/b, n/c) is the best nesting vector. Now the Green's function of an USDW is 
obtained as 

G _1 (k, iuj n ) = iuj n - £ k P3 ~ A(k)pi<7 3 , (9) 

while for UCDW 03 has to be replaced by one. We note at this point, that we assumed a real order parameter 
A(k), as in the absence of impurities and pinning the phase is unrestricted and therefore can be chosen to be zero for 
convenience. Here pi, (ai) are the Pauli matrices acting on momentum (spin) space respectively, while the linearized 
spectrum of the highly anisotropic electron system around the Fermi energy is £k = f ^—p = VF(k x — kF) — 2tb cos(bk y ) — 
2t c cos(ck z ). The order parameter is either independent of the momentum, which is the case of a conventional DW, 
or it can have four different type of wavevector dependence (A(k) = Acosbfcj,, A(k) = Asinbfcj,, A(k) = AcoscA: z , 
A(k) = Asincfc z ) as discussed in detail in Ref. |22|. Henceforth, without the loss of generality we can assume a k y 
dependent gap to be open, since the two perpendicular directions are equivalent in our model. 
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III. QUASIPARTICLE CONTRIBUTION 

Making use of the anisotropic nearest neighbor tight-binding band structure and Eq. © , the Raman tensor becomes 
diagonal with the same cosine functions being in the diagonal that appear in ek- Since the band structure belongs 
to the completely symmetric irreducible representation A g of the pointgroup of the lattice D 2 h, it follows that this 
is similarly true for every component of the vertex 7k- Our model is therefore only capable of describing the Raman 
spectra belonging to the A g symmetry channel, i.e. the spectra measured in x — x, y — y and z — z scattering 
geometries. In order to generalize the present treatment to incorporate the possibility of scattering effects with 
perpendicular polarizations, for instance the x — y geometry, nonvanishing offdiagonal components of the Raman 
tensor are needed. Particularly in the x — y geometry, the inclusion of a second nearest neighbor hopping term in the 
a — b plane in the one particle energy can account for finite absorption. It can be readily shown that in an orthorombic 
lattice the spectrum obtained in the x — y geometry belongs to the B\ g representation. We return to this point at 
the end of this section. 



A. Raman spectra with A g symmetry 

Coming back to the Green's function in Eq. J!|J), the quasiparticle contribution to the Raman susceptibility for 
(U)DW can be written as the sum of three terms corresponding to the three different polarization directions 

X 77 (q, iu n ) = -7377 Tr(r(k)G(k, iw m )T(k - q)G(k - q, iu m - iuj n )) 

PV ktt (10) 

X 77 ~~~ X"yy ~~~~ 'X.'yy J 

where T(k) is a four by four diagonal matrix with the elements (7k, 7k-Q> 7k, 7k— q) m the diagonal. For simplicity, 
we shall limit our analysis to q = (q x , 0, 0) (i.e. wave vector pointing in the quasi-one-dimensional direction). For the 
retarded correlation functions we get 



I 2 

\ u .,,^-^ : MOt^^^a fl-^^^-F^, (lib) 



X77 (£,c) = 2 5 (0) 7 2 \ M 2 7^J^(1 - 4A 2 F 2 ) + A 2 ( + co 2 F 2 - f ~ ^ 2 F 4 )}, (I In ! 



e 2 


-UJ 2 
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-LU 2 




e 
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T 2 



X 77 , cos (^) = 4g(0) 7 2 i 2 -^ 1 - ^j-F 4 ) , (11c) 



4 .2 A 2 

X UZ,u) = 4 ff (0) 7 2 t 2 -^ ( 1 - ——F 2 I , (lid) 



e-^ 2 \ $ 

where g(0), (p(0)) is the density of states at the Fermi energy in the normal state per spin per unit volume (per unit 
cell), £ = VFq x , /J- is the chemical potential and 7x = ma 2 (e x e s x ) , 7j/ = mb 2 (e y ey), 7 Z = mc 2 (e*ej). Furthermore P is 
the relevant coupling responsible for the DW formation whose detailed form can be found in Ref. [22], while 

t M 2^ 1 f°° f 2 * f@E\ N sir/%) , 

F n = (e - c 2 ) — / / tanh V 77 Re . Ky J^=dydE (12) 



2^ ./() ./() V 2 I D / r .> \ ■ .;,,2, 



'S 2 - A 2 sin%) 

TV = (£ 2 -w 2 ) 2 -4£ 2 (£ 2 +w 2 )+4A 2 £ 2 sin 2 (y) 
D = N 2 - QAE 2 uj 2 ( 2 {E 2 - A 2 sin 2 (y)) 

is the -F function that shows up in the correlation functions of conventional DWs with constant gap |3l| , as well as in 
the unconventional DWs [12 • Eqs. ifTT|) correspond to single bubble diagrams with self energy corrections due to the 
order parameter of the condensate. 

Since the momentum transfer of scattered light is small compared to the Fermi wave vector, we are only interested 
in the £ — ► long wavelength limit. Taking the imaginary part of the obtained susceptibilities according to Eq. 
at zero temperature we find for conventional DW 

Inw^.conv = 2 3 (0) 7 2 A 2 ^ReVx 2 - 1, (13a) 
Inw^.conv = M0H,A , e ^Re J_ , (13b) 
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FIG. 2: Raman spectrum of an UDW for y — y polarization with A(k) 
panel). Insets: the same spectra in a conventional DW. 




lo/2A 

Asin(&fc I/ ) (left panel), with A(k) 



Acos(6fc a ) (right 



while for unconventional DW we get 

2 5 (0)7^A 2 



Imx^ 7 = 



Im X 77 ,sin 



3.T 



2(x 2 - l)K(x) - (x 2 - 2)E{x), x<\, 
x(x 2 - l)K{l/x) - x(x 2 - 2)E(l/x), x>l, 



3x 

8g(0)7 y 2 ^ 
3x 



tH J (l-x 2 )K(x)-(l-2x 2 )E(x), 



x < 1, 



2x(l - x 2 )K{l/x) - x{l - 2x 2 )E(l/x), x>l, 



(x 2 + 2)K(x) - 2(x 2 + l)E(x), 



x < 1, 
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x(2x 2 + l)K{l/x) - 2x(x 2 + l)E(l/x), x>\, 

K{x)-E(x), x<l, 
x(K(l/x) - E(l/x)), x>l, 



(14a) 
(14b) 
(14c) 
(14d) 



where x — ui/2A and K(x), E(x) are the complete elliptic integrals of the first and second kind, respectively. The 
response functions for finite T are obtained simply by multiplying Eqs. l|13fl and (|14l) by tanh(w/4T). The spectra 
are shown in Figs. ^ an d0 

For conventional DW - due to the equivalence of the b and c crystal directions - not surprisingly we obtain the same 
response for the y — y and z — z polarizations showing the usual inverse squareroot divergence at 2A. It is worth 
mentioning, that for the chain polarization x — x the divergent peak is suppressed and transformed to a squareroot 
edge at the same threshold due to the vanishing vertex on the Fermi surface. It turns out that the vertex responsible 
for the scattering intensity is (7k — 7ic-q) 2 m our Nambu notation, which in this particular scattering geometry is 
proportional to that clearly vanishes on the Fermi surface. 

In contrast to conventional DW, in an UDW there are linenodes on the Fermi surface |22j . giving rise to arbitrarily 
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small energy nodal-excitations. It follows that the scattering intensity is finite for frequencies smaller than the 
maximum optical gap 2A. Furthermore, in the y — y geometry the interplay between the cosinusoidal vertex and 
the order parameter - either sinusoidal or cosinusoidal - results in two qualitatively different lineshapes (see Fig. [5J . 
The clear singularities in the spectra at 2 A for UDW that appear in Figs. ^ and [21 are of logarithmic type and are 
caused by Van Hove singularities in the quasiparticle density of states :22j. However this peak in Fig. |2 (left panel) 
is suppressed because of the zero Raman vertex at the gap maximum. 

The low frequency power law behavior is also characteristic for systems with point- or linenodes on the Fermi 
surface. In particular we have 

Imx*> - 0) = 2 5 (0) 7 2 A 2 ^ (^) 3 + 0(u; 5 ), (15a) 
In< 7 ,si> - 0) = 2 5 (0) 7 ^^ + 0(u; 3 ), (15b) 
Im^ 7)C( > - 0) = 2 5 (0) 7 ^ (^) 3 + 0(u% (15c) 

Im X ;> 0) = 2 5 (0) 7 ^!^ + 0(cj 3 ). (15d) 

All these important features of the Raman response makes the Raman experiment to be a relevant and fruitful probe 
in identifying the magnitude and symmetry of the energy gap. Similar analysis |33j contributed to the establishment 
of the ci-wave nature of the order parameter in HTSC. 



B. Raman spectra belonging to B\ g , B2 g and Bi g symmetries 

We have already pointed out in the beginning of this section that the choice of the nearest neighbor tight-binding 
band structure for the one particle energies is only sufficient to describe the Raman response in the A g symmetry 
channel. Now we extend the previous analysis with the inclusion of second nearest neighbor hopping terms in the 
a — b, a — c and b — c crystal planes, respectively. This extension on one hand can be considered as the simplest 
natural and physically motivated step towards the treatment of more realistic one particle energies, on the other hand 
it is sufficient to explain the Raman spectra measured with perpendicular polarizations that belong to the other three 
irreducible representation of the pointgroup. Namely, we add to the following extra term 

<5ek = At xy cos(ak x ) cos(bk y ) + 4t xz cos(ak x ) cos(ck z ) + 4t yz cos(bk y ) cos(ck z ). (16) 

The first term will clearly give nonzero offdiagonal component in the Raman tensor (^ xy ~ sm(ak x ) sm{bk y )). This 
function of the wavevector k is the most simple basis function belonging to the B\ g representation of £>2/i- Similarly 
the second and third terms of Eq. (|16fl yield nonvanishing contributions in the other two offdiagonal positions of the 
tensor, and the corresponding functions belong to the remaining two representations B 2g and B 3g , respectively. 

Now making use of the formalism and notations introduced in the previous subsection, the Raman spectra in the 
x — y scattering geometry - labelled by the B\ g symmetry - are obtained as 

{p% / i 4A 2 lj 4 \ 1 

at %-^ji (! - 8A " F ^ + ^ [tffip + <^ F * " {u) 2_ e) 2 F e) j . ( 17a ) 

X^cos&u) = ig^l Ltl v --^-(l-8A 2 (F 2 ~F i )) 



1 ->,„ 4A 



2 cj 4 



' 2pA2 ( W)p +uj2{F2 Fi) ~ {^-e? {Fi ~ Fe) ' ^ ' 1 171,1 



In the x — z geometry we have 



X^(^) = 4g(0) 7 L { a&^-^l - 4A 2 ^) + (3' A 2 + u?F 2 - ^ _ ^ 2)2 F 4 ) \ , (is, 



and finally in the y — z geometry we get 



Xr^infcw) = 8 9(°)7^^37j2 (1-8A 2 F 4 ), (19a) 
xftcciM = 8g(0)7, 2 z ^72^ (! " 8A2 (^ ~ F A )) . (19b) 
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Here F n is given by Eq. l|12f) . the superscripts sin and cos in Eqs. (|17(l and (|19f) denote the wavevector dependence of the 
order parameter A (k) as in Eqs. (|TT)) . Furthermore a = 4sin 2 (a/cF), (5 = (t xy /t a ) 2 cot 2 (ak F ), (3' — (t xz /t a ) 2 cot 2 (ak F ), 
-f xy = mab, = mac, and j yz = mbc. In the long wavelength limit £ — > and x Bza clearly vanishes, while x Blg 
and x B2a remain finite. The spectra emanating from the latter two are qualitativly not different from Imx^ 7 plotted 
on Fig. ^ left panel, therefore are not shown here. Nevertheless, the above considerations indicate, that assuming a 
more realistic band structure would not significantly alter the Raman lineshapes obtained in this section. 



IV. ELECTRON-ELECTRON INTERACTION, RPA SERIES 



The previous section dealt with the Raman response function in the one bubble approximation, i.e. the effect of 
interaction is taken into account in the self energy only. Now we turn our attention to vertex corrections at the 
RPA level, since the short and long wavelength components of the interaction may give rise to collective modes and 
Coulomb screening, respectively. 



A. Collective excitations in the Raman response 



Following Ref. [32] , the short wavelength component of the electron-electron interaction favoring a sinusoidal gap 
in the k y direction, namely A(k) = Asin(6fcj,), is given by 



N ~ 



5- a ,<7< {2J y sm(bk y ) sm(b(k' y — q y )) — 2F y sm(bk v ) sin(&fc^)) 
+ $<t,it'{Jv - V y ) sm{bk y ) sm.{b{k' y - q y )). 



(20) 



The detailed form of the whole interaction responsible for the density wave formation with the relevant couplings (P, 
see Eq. I|21|l) can be found in Ref. [22] • Here we shall continue with the assumption we made in Section II, namely 
that as we enter the low temperature phase (LTP), a gap varying in the k y direction opens up first and persists all the 
way down to zero temperature. Moreover we fix its functional form to be sinusoidal. All the calculations we present 
here can also be done with cosinusoidal gap without any relevant changes. 

In the small momentumtransfer limit (q — ► 0), using the spinor introduced in Eq. ©, the interaction operator 
corresponding to the matrix element in Eq. I|2U|) can be recast as 



H iat = -^*t( k + q)^(k)*(k)*t( k ' _ q)^(k')*(k'), 



(21) 



where i = c,s for UCDW and USDW, respectively. A c = p\ sin(6fc y ), A s = pia 3 sm(bk y ) and the detailed form of the 
couplings P c for unconventional charge- and P s for spin-density waves are given in Ref. |22j. With this, the correlator 
of the "effective density" in the RPA reads as 



X77 



XAij ' 



PiV c o 

I- — 7-Xa % a z Xa z1 , 



(22a) 
(22b) 



where V c is the cell volume and now the zero superscript denotes the one-bubble result. X77 1S already given in 
Eq. JTIJ. In addition we obtain 



Am 



(u/ - e)F2 - 4A^F 4 , 



■ lo 2 F 2 - 



oj 2 -e 



;F 4 



(23a) 
(23b) 



It is readily seen from Eq. I)23b|) that the contributions from the y—y and z—z scattering geometries do not appear in the 
nondiagonal XyA susceptibility, allowing only one Dyson series to develop, the one which dresses the chain polarized 
response, x — x. That is, taking into account the collective degrees of freedom of the DW, namely the fluctuation of the 
phase and amplitude of the order parameter around its mean field value, only this latter spectrum gets renormalized, 
while the former ones with polarizations aligned perpendicular to the quasi-one-dimensional direction retain their 
single-particle form. This situation is similar to the DC or optical conductivity of density waves, where the phase 
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FIG. 3: Left panel: the real (solid line) and imaginary (dashed line) parts of u> 2 F2 — 4A 2 _F4 at £ = and T = 0. Right panel: 
the RPA Raman spectrum of an UDW for the chain-polarized scattering geometry, x — x at T = 0. Inset: the same RPA 
spectrum in a conventional DW. 



mode of the condensate contributes only to the chain direction as well 32]. Solving the coupled RPA equations in 
Eq. , for the full susceptibility we have 

PiV c {x° Ai 



X77 = X° 7 + ^ V4 ■ (24) 



0,x 

A-yy 



From this the Raman intensity is obtained as 

Im X * 7 = Im- {i . 

4 XA,A, (25) 

where A = g(0)PiV c is the dimensionless coupling. The plot is shown in Fig. at T = 0. The nature of the collective 
mode can be simply explored as usual by looking at the poles of the response function in Eq. (|24[l . With Eq. I|23a|) 
the task is reduced to finding the roots of 

(lu 2 - £ 2 )F 2 - 4A 2 F± = 0. (26) 

We note here that a similar expression ((w 2 — C 2 )^2 = 0) appears when considering the density-density correlator 
of an (U)DW [32^. from which one obtains - for both conventional and unconventional DWs - the well known phason 
dispersion uJ 2 = £ 2 [34J- In our case for a conventional DW, since F2 — F4 (in Eq. (|12l) sin(y) has to be replaced 
by 1, because neither the interaction nor the order par ameter depends on the wavenumber) we reobtain the gapped 
dispersion of the amplitude excitation uj 2 = 4A 2 +£ 2 j34j, which is known to couple to the Raman experiment p] l35l36j . 

Now turning our attention to the unconventional situation and remembering the result of the previous case, we 
look for the root on the real frequency axis around 2 A. Therefore - after analytic continuation of the F n 's for real ui - 
we plot the lhs. of Eq. l|^rj|l versus frequency in the £ — > limit, see Fig. left panel. It is clear that - unlike what is 
found in Ref. [37j for a c£-density-wave - there is no zero in either the real or the imaginary part at 2 A, or at any other 
real frequency. There is no indication, that there would be a root of Eq. 126(1 for complex frequency either. It follows 
that in contrast to conventional systems, in UDWs although the Raman intensity shows considerable renormalization 
due to electron-electron interaction with respect to the one-particle form, there is no clear, particle-like mode with 
infinite lifetime with which we could identify the peak around w = 2 A, see Fig. |3 right panel. We can say however 
that the Raman vertex couples to the amplitude mode of the condensate, overdamped because of the low energy 
excitations. 
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B. Coulomb screening 

In the case when the light produces charge fluctuation in the electron gas, the coupling to the long-range Coulomb 
forces reduces the scattering rate. Therefore it is useful to treat these forces separately. In the usual RPA approach 
the screened Raman susceptibility reads as 

X 77 = X77 - X71 (Y - VxnV + . . .) xi 7 = X77 h 2 — Xu, (27) 

X11 X11 

where V = Aire 2 /q 2 , X77 is already calculated in Eq. (fTTfl . xW = Xn/(1 + Vxn) is the screened density correlator 
with xn being the single-particle contribution in unconventional charge- and spin-density waves given in Ref. |32| . 
The pole of the screened density correlator leads to the plasmon mode [38|. which however does not affect the low 
energy behavior. Furthermore, the nondiagonal X7I term is expressed as 

X71 = -0(O)7«(2A* + u) - 4A 2 F 2 ). (28) 

At this point it is important to call the attention to the fact, that similarly as in the previous subsection dealing 
with the short wavelength component of the interaction, here again we see that only one component of the whole 
Raman vertex survives in Eq. 1|28|) . namely the x — x component. It means that in principle only the Raman spectrum 
with the incoming and scattered polarizations aligned in the chain direction can be screened by the Coulomb forces. 
Nevertheless, since both the xu ancl X71 correlators are quadratic in momentum in the long wavelength limit, thus 
the coupling between these quantities is not that strong to modify the zeroth order term x 77 in Eq. 1|27|) . The above 
calculations therefore lead us to the conclusion that the Raman response is not affected even if we take into account 
Coulomb screening in RPA. This is due to the vanishing average of the Raman vertex on the Fermi surface. 

V. CONCLUSIONS 

We have investigated theoretically the electronic Raman scattering in quasi-one dimensional interacting electron 
systems with density wave ground state. Mean field treatment of conventional as well as unconventional density waves 
in pure systems has been applied in order to determine the Raman intensity in various scattering geometries. We 
have found distinct, characteristic lineshapes especially in the unconventional situation, depending on the particular 
momentum dependence of the density wave order parameter. We conclude, that the Raman experiment could serve 
as a valuable tool in identifying materials supporting unconventional density waves, and in specifying their particular 
gap structure. We have also considered Coulomb screening, and we found it ineffective due to the negligible coupling 
of density fluctuations to the Raman vertex in our nearest neighbor tight-binding model. Collective contributions to 
the Raman response appear only in the x — x scattering geometry (light polarized in the chain direction) . This is due 
to coupling to the amplitude mode of the condensate. This mode is overdamped in the unconventional case, since 
decay to low energy excitations is possible. 
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